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Abstract 



£> ' The QCD sum rules for spin-dependent nucleon-nucleon interactions are for- 

mulated and their physical implications are studied. The basic object of the 

' study is the correlation function of the nucleon interpolating field, where the 

00 ' 

matrix element is taken with respect to the one-nucleon state. By means of 
the dispersion relation, the correlation function in the deep Euclidean region, 
where it is expressed in terms of the nucleon matrix elements of the quark- 



gluon composite operators by using the operator product expansion (OPE), 



^ ' is related with its integral over the physical region. The dispersion integral 

S-T 



of the correlation function around the nucleon threshold is investigated in 
detail. It turns out that the integral can be identified as a measure of the 
nucleon-nucleon interaction strength, which is proportional to the scattering 
length in the small scattering length limit and to one half of the effective 
range in the large scattering length limit. New operators, such as q^^ 75(7, 
<j75°~fivQ, must be taken into account in the OPE of the correlation function. 
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Therebehavior operators do not vanish when the matrix element is taken 
with respect to the spin-nonaveraged one-nucleon state. The Wilson coeffi- 
cients of such operators are calculated. The sum rules obtained in this man- 
ner relate the spin-dependent nucleon-nucleon interaction strengths with the 
spin-dependent nucleon matrix elements of the quark-gluon composite oper- 
ators. The sum rules imply that the interaction is stronger in the spin-triplet 
channel than in the spin-singlet channel, but that the spin-dependence of the 
nucleon-nucleon interactions is rather small. In the spin-singlet channel the 
calculated strength is in qualitative agreement with the empirical strength, 
which is estimated by the empirical low energy scattering observables. 

I. INTRODUCTION 

One of the most important goals of investigating strong interaction physics 
is to understand the of hadrons and hadronic interactions on the basis of 
quantum chromo dynamics (QCD). Shifman, Vainshtein and Zakharov pro- 
posed the method of the QCD sum rule, which provides us with a framework 
to investigate the properties of hadrons in a model-independent way. || This 
method has been successfully applied to the study of the masses, decay con- 
stants, magnetic moments and other properties of various hadrons. Q 

Recently, the present authors extended the QCD sum rule to the investi- 
gation of hadronic interactions. Q In Ref. 3) a nucleon-nucleon system was 
studied as a typical case. The correlation function of the nucleon interpolat- 
ing field, whose matrix element is taken with respect to the one-nucleon state 
and averaged over the nucleon spin, was considered. It was noted that the 
correlation function has a second-order pole at the nucleon on-shell energy 
as a function of the energy associated with the interpolating field and that 
its coefficient is the T-matrix for the nucleon-nucleon (NN) scattering. As- 
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suming that the dispersion integral is dominated by the pole term, sum rules 
were derived which relate the spin-averaged NN scattering lengths with the 
spin-averaged matrix elements of the quark-gluon operators with respect to 
the one-nucleon state. The obtained NN scattering lengths are of the order of 
several fm, which is rather large in the strong interaction scale, but are rather 
smaller than the experimental values. The formalism was further applied to 
other hadron-nucleon systems. j|] || 

The following point, however, remained unclarified in those works. In 
the analysis of the QCD sum rule for the hadron in the vacuum, the cor- 
relation function of the hadron interpolating field, whose matrix element is 
taken with respect to the vacuum, is considered. The imaginary part of the 
correlation function consists of a pole term corresponding to the ground state 
and a continuum term corresponding to the excited states. Under the Borel 
transformed dispersion integral, the continuum contribution is exponentially 
suppressed compared to the ground state contribution due to the energy differ- 
ence between the ground state and the continuum threshold. For this reason 
the sum rule analysis is expected to be insensitive to the detailed form of the 
continuum, so that the continuum is usually parametrized in a very simple 
form. When one deals with the hadron correlation function, whose matrix 
element is taken with respect to the one-nucleon state, the situation is dif- 
ferent. The energy of the continuum threshold is not higher than the pole 
energy. Therefore, it is not clear if the Borel transformed dispersion integral 
is really dominated by the pole term or not. 

A related question is the following. It is known that there is a loosely 
bound state in the spin-triplet nucleon-nucleon channel and an almost bound 
state in the spin-singlet nucleon-nucleon channel. If there is a zero-energy 
bound state, the scattering length diverges. Therefore, the NN scattering 
lengths are expected to be very sensitive to the NN interaction strength. 
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On the other hand, it is hard to believe that the nucleon matrix elements 
of the quark-gluon operators would be very sensitive to the NN interaction 
strength. It seems strange that the sum rule relates these two quantities of 
very different natures. 

Another point is that only the spin-averaged sum rules are obtained from 
the spin-averaged correlation function. The NN channel is special in the 
sense that selecting the isospin channel automatically selects the spin state. 
In Ref . [|| , the sum rules for the spin-triplet and singlet scattering lengths are 
obtained by combining isospin states. As far as the pole term is concerned, 
the above selection rule is correct, but it does not hold for the continuum. 
The sum rule in Ref. 3) is valid only if the pole term is dominant. Therefore, 
it is more desirable to construct the spin-dependent sum rules. 

In this paper we consider the spin-dependent correlation function of the 
nucleon interpolating field, where the matrix element is taken with respect 
to the spin-nonaveraged one-nucleon state. The purpose of this paper is two- 
fold. First, we extend the procedure of the sum rule to the case of the spin- 
nonaveraged correlation function. Second, we show that the dispersion inte- 
gral of the correlation function around the nucleon threshold can be regarded 
as a measure of the nucleon-nucleon interaction strength. As a result, we de- 
rive sum rules which relate the spin-dependent AW interaction strengths with 
the spin-dependent nucleon matrix elements of the quark-gluon operators. 

II. FORMULATION 

A. Physical content of the correlation function and Borel sum rules 

Consider the spin-dependent correlation function, IL(qps), 
U(qps) =~ij d 4 xe iqx (ps\T(iP{x)ip(0))\ps}, 
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where \ps) is the one-nucleon state with momentum p and spin s (p 2 = M 2 , 
s 2 = — 1 and ps = 0, where M is the nucleon mass) normalized as (ps\p's') = 
(2-7r) 3 (5 3 (p— p')5 ss i and ip is the normalized nucleon field operator (0\ip(0)\ps) = 
u(ps), where u(ps) is a positive energy solution of the free Dirac equation 
for the nucleon. In this paper, momentum with " represents the nucleon 
on-shell momentum. Later, the normalized nucleon field, ip, is replaced by 
the unnormalized nucleon interpolating field (quark-gluon composite field), 
r). The following discussion, however, holds as it is for the interpolating 
field, except for the normalization. Naively, the dispersion relation for the 
correlation function, Tl(qps), is written as 

1 f°° , 1 , 
U(qps) = / dq r ——lmU(qps), 

vr J-oo qo ~ % + it] 

where q' = (q'o,q)- Throughout this paper, whenever we take the imaginary 
part of a quantity, we approach the real energy axis from above in the complex 
energy plane. Therefore, strictly speaking, Imll is the imaginary part of the 
retarded correlation function. The QCD sum rules are obtained by evaluating 
the left-hand side of Eq. ([!]) by the operator product expansion (OPE) and 
expressing the right-hand side in terms of physical quantities. 

Let us consider the singularities of H(qps) as functions of qo. In the 
complex qo plane, U(qps) has a branch cut from the lowest B = 2 continuum 
threshold to the right and another branch cut starting from the lowest B = 
continuum threshold to the left. In addition, U(qps) has second-order poles at 
qo = ± y q 2 + M 2 = ±Eq whose coefficients are the iV\/V and NN T-matrices 
T + and T_, respectively: 

T + {qrps; qrps) 

= -i [ d 4 xe iqx J^u{qr)(4- M){ps\T{^{x)$(Q))\ps){i- M) 
J y Eq 

= (qo - E q ) 2 ^-u(qr)U(qps)u(qr), 
Eq 




T-(qrps; qrps) 



f / A/f / A/f 

i / d A xe iqx J— v(qf)(4-M)(ps\T(ip(x)$(0))\ps)(4-M)J— v(qf) 
J V -Eg V E a 



M 



M 
Q 



(q Q + E q ) —v{qr)U(qps)v(qr). 



Here q = (qo, —q), f = (vq, —r) and v(qf) is the negative energy solution of 
the free Dirac equation for the nucleon. 

In order to take out the pole contribution from Imll(qps) it is convenient 
to define off-shell NN and NN T-matrices by 



T + (q'r'p' V; qrps) 
—i \ d A xe 



M){p's'\T^(x)^))\ps){i ~ M)J^-u(qr), 

an 



T_(q'r'p's'; qrps) 

-i J d*xe**J^v(qr)U- M)ip's'\T^{x)m)\m^ -M)0-v^f'). 



Note that Eq. (|2|) is just a definition of the T-matrix off the mass shell, but 
the LSZ reduction formula shows rigorously that it is the T-matrix on the 
mass shell. 

In order to separate the contribution from the poles at qo = Eq and 
qo = —Eq, we introduce the projection operators A + and A_ by 

A+(qs) = u(qs)u(qs) = t^L— tp^, 
A-{qs) = v(qs)v{qs) = , 

which have the properties 

A 2 + (qs) = A+(qs), 
A 2 _(qs) = A_(g a ), 

A + (qs)A + {qs) = A_(gs)A_(gs) = A + (ga)A_(gs / ) = 0. 
Then we define the projected correlation functions by 
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M 

U + (qrps) = — tv {A + (qr)n(qps)} , 
&q 
M 

U-(qrps) = ti {A- (qr)Il(qps)} . 

Eg 

The projected correlation functions are related to the off-shell T-matrices 



as 



n±( ^ s) = (^F' 

where T±(qrps) = T±(qrps;qrps). Clearly, TL±(qrps) has a second-order pole 
at go = ±Eq but not at go = TEq. 

Naively, the dispersion relation for the projected correlation function, n_|_, 
is given by 

1 f°° , 1 , 

U + (qrps) = / dq'n ; Imll+(q rps). 

vr J-oo qo-% + iil 

Formally, the imaginary part of the correlation function is written as 

lmll + (qrps) = Im ^ReT + (qrps) + Re ^lmT+(qrps) 

(qo ~ E q + irj) (q Q - Eg + irj) 

Pf 

= ir5'(qo — Eq)ReT + (qrps) H ^ m, E + (qrps) 

(q - Eg) 

i Pf 
= w6 {q - E q )t - ir5(qo - E q )u + -^lmT + (qrps), 

(qo-E q ) 

where 

t = ReT + (qrps)\ qo=Eq , 
d 

u = — — ReT + (qrps) 

d( 10 qo =Eg 

However, as it will turn out, when q = and p = 0, the integral of the second 

term is divergent because u is divergent, and the integral of the third term is 

3 

also divergent because it behaves as (go — M)~z in the vicinity of go = M. 
Therefore, Eq. (@) is ill-defined. 

Instead of Eq. (||), we consider the dispersion relation for ^ ^ R(qrps), 

4 ~ M„, „ . 1 f°° , , 1 r (4' ~ M^, . . , 

U(qrps) = / dq Q ——Im \ — II(g rps) 

qo vr J-oo qo-q + *V I <?o 



or for q ° ^ Yl^(qrps) in terms of the projected correlation function, 



— ^-U + (qrps) =~ [ dq \ , Im j % Eq Ii + (q'rps) \ . 

Qo 7T J_oo q - q' + irj { % J 

Symmetrizing Eq. (||) we obtain 
go __g n ( p s j + ( _ q ^ = - f dq' - — ^ aiq'o ~ E q )lmU + (qrps), 

where 

P 

(qo - Eg)lmll + (qrps) = -7rd(q - E q )t -\ — ImT + (qrps). 

qo - Eq 

Now Eq. (||) is well-defined when q = and p = because u does not 
appear and the second term behaves as (qo — M)~2 in the vicinity of go = M. 
Applying the Borel transformation, 

(ql) n ( d 



L B = lim . , 



to both sides of Eq. (||), we obtain 
-qo — Eg 



L B 



2qo 

Loo 6XP ( ~ M 1 ) ( q '° ~ E Q) lmU +( < l' r P s )' 



U + (qrps) + (q -> -q ) 



7T 



•y'2 
?0 

L B \ ±vI B, 

where Mb is the Borel mass. In order to derive the Borel sum rules we must 
evaluate the left-hand side by the OPE and parametrize the right-hand side 
in terms of physical quantities. Now, the question is how to parametrize the 
integrand of the right-hand side of Eq. (|7|) . 

Let us recall the QCD sum rule for the nucleon in the vacuum, where the 
imaginary part of the correlation function has the form 

Imn + (<?) oc -7n5(go - Eg) + a(q). 
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The first term is the contribution from the nucleon pole term, and the second 
term is due to the excited states. Borel transformation on the dispersion 
integral of Imil_|_ (q) gives 

r 1 



L B 



— IL + (q) + {q 
2<?o 



1 1 dQ'o^e-^/ M hmU + (q') 



where the second term, the contribution from the excited states, starts at the 
continuum threshold, to (u > Eq), and is exponentially suppressed compared 
to the first term. For this reason, it is possible to use the rough model of the 
hadron continuum, 

Imn+(g) = -A 2 7T% - E q ) + {8(q - cu ) + 9(-cu - q )} Imn° PE (g), 

where H+ PE is the asymptotic form of the correlation function in the OPE, 
ujq is the effective continuum threshold, and the normalization constant A is 
explicitly included ({0\r](0)\ps) = Xu(ps)). 

Let us turn to the problem at hand. As an extension of Eq. ( |l~7| ) one might 
parametrize (go — Eq)lmll + as 

{q - E q )ImU + 

= -A 2 7T% - E q )t + {6(-q - co.) + % - ^+)} (go " E q )Imn° PE , 
by approximating the second term of the right-hand side of Eq. (0) by its 



asymptotic form. However, the second term starts at qo = co = y 4M 2 + q 2 ■ 



M (u! < \J M 2 + q 2 ), and it is not exponentially suppressed compared to the 
first term. Therefore, one cannot justify Eq. (||). One has to know the behav- 
ior of the second term around the threshold. For this purpose it is important 
to note that the off-shell optical theorem holds for T. When the center-of-mass 
energy is above the threshold of the NN channel and below the threshold of 
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the next channel, only the NN states contribute in the intermediate states, 
and the off-shell optical theorem is simplified as 

^y^^^(2vr) 3 (5 4 (p + q - p n - q n )T + (qp; q n p n )T + (q nPn ; qp). 



In order to simplify the notation we introduce the scattering amplitude / by 

/ / ^ ;1 /V /2 / / 

f(qp; qp) = — T + (q p ; qp), 

where a = 9u f" and u! = . Moreover, we go to the center-of-mass frame 
(q + p = q' + p' = 0) and restrict ourselves to the s-wave. We define three 
scattering amplitudes, /o, f\ and /2 as 

fo(k) = f(q'p';qp), 

where \p\ = \q\ = \p'\ = \q'\ = k, p Q = q = p' Q = q' = \/M 2 + k 2 , 

fi(k)=f(q'p';qp), 

where \p\ = \q\ = k, \p'\ = \q'\ = 0, po = qo = V M 2 + k 2 , p' = M, 
q' = 2\/M 2 + k 2 - M, and 

f2(k) = f(q'p';qp), 



where \p\ = \q\ = \p'\ = \q'\ = 0, p = p' = M, q = q' = 2^M 2 + k 2 - M. 
It is well known that the on-shell scattering amplitude /o has the form 

/o(fc) = -ik + k cot 5 = + I + IrA; 2 + 0(fc 4 ) ' 

where a is the scattering length and r the effective range. Similarly, the 
off-shell scattering amplitude /2 has the form 

f2(k) = ' 



i{—k + bk 3 + 0(k 5 )} + {i + \fk 2 + 0(fe 4 )} ' 
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which can be shown as follows. First, since the discontinuity of the T-matrix 
along the real energy axis is proportional to the imaginary part of the T- 
matrix, the Taylor expansion of the real and imaginary parts of the scattering 
amplitude includes only even and odd powers of k, respectively. Second, /o, 
fi and f2 coincide on the mass-shell (k = 0), /o(0) = /i(0) = /2(0) = a. 
Therefore, we have 

Re T7lA =- + \rk 2 + 0{k A ), 
fo{k) a 2 

J2{k) a 2 



Third, from Eq. ( |1Q[) the following relations hold, 

lmf (k) = k\f (k)\ 2 , 
lmf 2 (k) = k\f 1 (k)\ 2 . 

Therefore, we have 

J_ Im/ (fc) _ 
fo(k) \fo(k)\ 2 

lm h(k)- \Hk)\ 2 ~ \Hk)\ 2 ~ + + ( } ' 

Equation (|Tl| ) follows from Eqs. (|l2] ) and (^). 

It should be noted that r is different from the effective range r, but f 
coincides with r in the limit a — > oo: 

f = r + O ^— ^ . 

This is shown as follows. Equations ( Jl2| ) and (13) indicate that both O(l) and 
O(fc) terms of l//o(fc) and l//2(fc), which are real and imaginary respectively, 
coincide with each other. Therefore, we have 

£S9-i + (*»). 

7o(«) 

Equation ([14]) is independent of a and therefore holds also in the limit a — > oo 
due to the continuity. Since l/fo(k) and 1/ f%(k) do not have the O(l) terms 



in this limit, they must coincide with each other up to the 0(k 2 ) terms, i.e. 
f = r + 0(±). 

From Eq. (Ill]) we have 



where k 



a + a 2 n + 0(n 2 ) qo < M 

1 + £) fl3fc2 + °( k4 ) 10 > M 

q < M 

a 2 k + 0(k 3 ) q > M 

-ik. One sees that 



Re/ 2 (k) 



Imf 2 (k) 



dqo 



ReT 4 



1 d (q + M 



q =M 



2ir dqo [ qoM 



Re/ 2 



oo, 



q =M 



and 



ImT_L 



1 qo+M 



Im/ 2 oc (q - M) 2 , 



2vr q M 

which make the naive dispersion relation, Eq. (|J), ill-defined. 

Having understood the structure of ImT + , we proceed to the integral of 
the right-hand side of Eq. ([/]), I, in the vicinity of qo = M, 

/ =4L^4 exp (^) w_M)Imn+ ' 

The integral, I, can be decomposed as 



I = I t + I c {+h 



In Eq. (15), the first term, It, is the threshold contribution, given by 



M 2 \ 4-n-a 



M 2 B 



exp 



Ml) M 



The second term, I c , is the continuum contribution, given by 



\ f dq'o-^exp (- 



J2 



Ml q' -M 



q' M J2 ' 



where 
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A: 



-%o-M). 



The last term, It,, is the bound-state contribution, which has to be taken into 
account if there is a bound state, given by 



Ml q' — M 



-2vr- 



q' M 



-Im/, 



poie 
2 



where 



Im/f' e = 



a / 1 



— ITT 



I dk V/2 
= itc5(k — Kq), 



k=iKQ 



6(k - kq) 



and iAto is the pole momentum, l//2(i«o) = 0. 

By performing the integral, the continuum contribution becomes 



-lI dq '°M- 



,/2 



exp 



Ml 



exp - 



B 

M^ 

Ml, 



-2vr 



q'o + M 



Ml) { q'o ~ M V ; ^ + + f + fc2 



fc 



47r|a| 



^ + m* (i + s + £) 

which is simplified in two limits of a as 



i exD 1-^-1 . 47r .. q 2 

eX P I M%) VT+M^b ' 



(a^O) 

m^ x PI-S)^I«K 1 -^)+°(h)- (a -oo) 
Similarly, the bound-state contribution becomes 

12 \ ( p 



1 2k 



^l 6XP l MlJU -M 



i 



/2 



1 



exp 



Ml) J -M 



2ir- 



q' M 
M + J 



Muj- 



—ire, 



where u/ = 2\/M 2 — Kq — M. In the limit a — > oo, 



K 



and the bound-state contribution is simplified as 
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h ^ ik exp (-^) S a i 1 - h) + ° G) • (a 00) 

Let us suppose that one can freely change the interaction strength of two 
nucleons and examine how the integral I should change as a function the 
interaction strength. When the interaction is weak, the scattering length is 
also small, and the integral / is dominated by It- 

I = h + h 

1 / M 2 \ 4vr ^ 2 

= T7T ex P 7TT TTfl + Old ). 

M B \ M B/ M 

As the interaction becomes stronger, the scattering length increases and the 
integral, I, also increases. As the interaction strength increases further, the 
scattering length eventually diverges when the bound state is just formed. 
Just before the bound state is formed, the integral I becomes 

I = It + Ic 

1 / M 2 \ f 4™ 4vra / r\ „ ( 1 \ 1 
1 / M 2 \ 2vrr /1\ 

and just after the bound state is formed, it becomes 



I = I t + I c + I b 




This shows that before and after the bound state is formed the integral is 
continuous, though the scattering length diverges with opposite signs. This 
observation leads us to conjecture that the integral around the threshold is a 
measure of the NN interaction strength. Based on this conjecture we define 
the NN interaction strength, a, by 
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1 / M 2 \ 4vra 
= "M| eXP \~Ml) ~W 

In the dispersion integral, the imaginary part of the correlation function, 
ImII + , contains the contribution from all possible intermediate states such as 
those of the AW, NNtt channels and so on. However, only the AW channel 
contributes around the threshold. We assume that the contribution from the 



AW state is taken into account by the form of the right-hand side of Eq. (16) 
and that the rest is approximated by the asymptotic form of the correlation 
function starting from an (effective) threshold, u+, for the B = 2 channels 
other than the AW channel and LO- for the B = channels: 

(go - M)Imn+ 

A 2 7t% -M)—a + {9(-q + % - ^+)} (go - M)lmR% PE \ 

A* 

where the normalization constant A is explicitly included. This is possible now 
because the contribution from states other than those of the AW channel is 
exponentially suppressed compared to the AW contribution. 



B. OPE of the correlation function and results 

Let us turn to the OPE. We take the interpolating field of the neutron 
as § 

ri(x) = e abc (d Ta (x)C lfl d b (x)) 75 7^ c (^), 

where C denotes the charge conjugation operator and a, b and c are color 
indices. We take into account all the operators of dimension less than or 
equal to four. We also include four-quark operators, they are of dimension 
six. The operators which involve the quark mass are ignored. In the OPE, 
the neutron correlation function, where the matrix elements are taken with 
respect to the one-nucleon state, is given as 
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n 



OPE 



(qps) 



4vr 4 



1 i 

7 ^[g 2 ln(-g 2 )vr 2 { - ~(d^d) N - -{u^u) N } 

2 2 
+<^<f ln(-g 2 )vr 2 { - -{d^ v d) N - -{vrf v u) N ^ 



1 lot 



O \ 7T p 



+ 



1 
1 

4^ 



+<z"ln(-<zV{ 

1 f 8 — - ~i i 
+<?^— 7r 4 |-(dd(i(i) w | 

q 2 ln(-q 2 )ir 2 {-(uu) N } + q^— — (uud<y^d) N } 



7 M 75 



r 5 1 i 

t ln(-q 2 )n 2 {-{d'y IJt 'y 5 d) N - -(u^^u) n \ 

2 2 
+g M <f ln(-g 2 )7r 2 { - -(dj v -y 5 d) N - -{u^ 5 u) N } 

+g ly ln(-g 2 )7r 2 | - ^{S[dj^ 5 iD u d]) N + ^(S[uj^ 5 iD u u\) N } 



+Q U ~2 71-4 { ~ y(^75^^^)iv} 



11 2 

+ ^75^[<? 2 ln(-<? 2 )7r 2 j - -(u7 5 f7^u}jv} + ^<fln(-g 2 )7r 2 { - -(«75^i/p«)iv}} 



In Eq. (18), D is the covariant derivative, S[A^B V ] = A^B U ^ — (traces) , 
where / i represents symmetrization over the Lorentz indices and —(traces) 
represents the subtraction of the trace terms. Here, (O)jv is the connected 
part of nucleon matrix element of O, (0) N = (N\0\N) - (N\N)(O) , where 
(0)o represents the vacuum expectation value of 0, and \N) = \ps). 

Substituting the projection of Eq. ( |i~8| ) into the left-hand side and Eq. ( [lT[ ) 
into the right-hand side of Eq. (j^), respectively, and splitting a into spin- 
independent and spin-dependent parts, a = a indep + a dep , we obtain sum rules 
for the spin independent and dependent interaction strengths, respectively, as 

M 2 \ 



\ 2 ^ n indep 1 exn 



M 2 B , 



4 ^4 { (C 2 MM B - C 3 M 2 ) [vr 2 { - 3(dU) N - (u^u) N } - ir 2 {uu) N 
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+ (C X M - C 2 M B ) [vr 2 {- l—G a ^G, 



'aft 



+ ■ 



16 



JV 



N 



+ 



3 

M 



4 i 

i(S[d7o-D d]> at + -i(S[«7oAH)Ar j 

;7r 4 (dddd) at + —^{uud^d) n \, 



and 



-A —a 

M 



dep 



Ml 



exp 



Mj 



l —r k { (C 2 MM B - C 3 M 2 B 



4vr 4 



x 7r 2 |-(d7 fc 7 5 d)jv - ^(n77 5 u)Ar} - ^7r 2 (u7 5 o- fc0 u)Ar 
+ (CiM - C7 2 M B ) 



7T 



M 



Ml 



—n 4 (ddd'y 5 ia k0 d) N - —TT 4 {uud-f k ^d) N j 



where 
1 



1 



exp 



Ml) 

( ojI 



+ exp 



M B 
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and $ is the error function. 

From the QCD sum rule for the nucleon in the vacuum in Ref. ||, the 
normalization constant A 2 is related with the vacuum expectation values of 
the operators according to 
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1 
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, 1 ^0 

M 2 B 2 Ml 



exp 
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The sum rule for the spin-independent part of the interaction strength, 
a mdep , is nothing but the sum rule for the NN scattering length in Ref. ||, 
which was obtained starting from the spin-averaged correlation function. The 
difference, however, is that the scattering length in Ref. || is replaced by the 
interaction strength in Eq. (|19[) . 

Equation ( |20|) is the new sum rule for the spin-dependent part of the inter- 
action strength, a p . The new sum rule provides us with a relation between 
the spin-dependent NN interaction strength and the nucleon matrix elements 
of the spin-dependent quark-gluon operators. Therefore, if one knows the nu- 
cleon matrix elements of the spin-dependent quark-gluon operators, one can 
predict the spin-dependent NN interaction strength, and vice versa. In this 
paper we investigate the first possibility. 

We now discuss the nucleon matrix elements of the quark-gluon operators. 
Dimension-three operators are quark operators, qq, q^y^q, (77^ 75(7 and cpi^o^q. 
The nucleon matrix elements, (qq)N and (qj^q)N, are spin-independent and 
have already been discussed in Ref. 7), while ((77^75(2) jv and (5750 ^q) n are 
spin-dependent and are written in terms of the axial charge, Aq, and the 
tensor charge, Sq, as 

The axial charge and the tensor charge are related to the structure func- 
tions H as 

Aq = [ dx[gx{x) + gi(x)}, 

Sq = dx[hi(x) — hi(x)], 
Jo 

where gi is the longitudinal quark-spin distribution in the nucleon, g\ is the 
antiquark-spin distribution, h\ is the quark-transversity distribution, and h\ 
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is the antiquark-transversity distribution, hi, together with the unpolarized 
quark distributions fi and gi, forms a complete set of twist- 2 structure func- 
tions. 

Recently Aq has received much attention in connection with the spin 
content of the nucleon. While a naive quark model picture suggests that the 
nucleon spin is carried by quarks, i.e. An + Ad = 1, it was experimentally 
found that Au + Ad << 1. || In Ref. 10), Aq for the proton was determined 
from recent EMC/ SMC and SLAC data with the SU(3) symmetry and the 
hyperon (5 decay as 

An = 0.83 ±0.03, 
Ad = -0.43 ± 0.03, 

at the renormalization scale Q 2 = 10 GeV 2 . There is also a lattice QCD 
calculation [11] with the result 



An = 0.638 ± 0.054, 
Ad = -0.347 ± 0.046 

at Q 2 = 2 GeV 2 . 

On the other hand, up to now no experimental information is available 
for 5q, since hi and hi cannot be measured by deep inelastic scattering. 
These quantities can, however, be measured by Drell-Yang processes, which 
are planned in a future RHIC experiment. However, 5q for the proton was 



calculated on the lattice [12| with the result 



5u = 0.839 ± 0.060, 
Sd = -0.231 ± 0.055 

at Q 2 = 2 GeV 2 . In this paper we use the lattice results both for Aq and 5q 
and ignore the Q 2 evolution of these matrix elements between Q 2 = 1 GeV 2 
and 2 GeV 2 . 
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Dimension-four operators are gluon operators, ^-G^G^, ^-S[G^G pi ,], 
and quark operators, S[qj^iD u q], qS{^^iD u )^q. The nucleon matrix ele- 
ments of the gluon operators are spin-independent and have already been 
discussed in Refs. 7) and 13). The matrix element of the quark operator, 
(S[q~'y ll iD v q\)N, is also spin-independent and is related to the unpolarized 



quark distribution /i, |14| while {qS{^^iD v )^^q) n is spin-dependent and 



By neglecting the operator including the quark mass, a\ can be related to the 
first moment of the longitudinal quark-spin distribution, g2, and it is given at 



tree- level by [15| 



a\ = — 2 / dxxg2{x). 
Jo 

Very recently, measurements of Q2 for the proton and the neutron have 
begun. |l6| |L7]] We have calculated a\ for the proton using the data for #2 
in Ref. 16) over the range 0.075 < x < 0.8 and 1.3 < Q 2 < 10((GeV/c) 2 
for the proton and those in Ref. 17) over the range 0.06 < x < 0.70 and 
1.0 < Q 2 < 17.0 ((GeV/c) 2 for the neutron with the results 

a? = 0.05 ± 0.04, af = -0.08 ±0.13. 

In addition to the above dimension-three and dimension-four operators 
we take into account the dimension-six four-quark operators, qqqq, qqq'J^.q, 
QQQJ[tl5Q an d qqqibO'fii/q, since four-quark operators are known to give the 
largest contribution among higher order operators in the QCD sum rule for 
the nucleon in the vacuum. M pi In the vacuum, the matrix elements of the 



four-quark operators are estimated by the factorization hypothesis; Q || i.e. 
it is assumed that the vacuum contribution dominates in the intermediate 
states: (0102)0 ~ (Oi)o(02)o- Similarly, for the nucleon matrix element, 
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we assume that the contribution from the one-nucleon state dominates in the 
intermediate states: 

(O x 2 ) N = {N\Oi0 2 \N) - <JV|J\O(0i0 2 >o 
„ (N\Oi\N)(N\Q 2 \N) 
(N\N) {°Vo{02)o{N\N) 

= (Ol)Ar(O 2 )0 + (Ol)0(O2>7V. 
Thus we assume {qqqq)N = 2(qq) {qq)N, (qqqj^N = (<?<7>o(<?7m<?>7V 5 

{qqqiiJ,im)N = {qq)o{qivnsq)N and {qqcftsa^N = (qq)o(qi5<y^q)N- 

For completeness we list here the values which were used in the calculation. 
The spin-independent proton matrix elements of the operators are 

{v)u) p = 2, (Sd) p = 1, (uu) p = 3.46, (dd) p = 2.96, 
i{S[u-i^D u u]) p = 222 MeV, %{S[d'y ll D v ^) p = 95 MeV, 
(^±G^Gn P = "738 MeV, (^S[G^°]) P = -50 MeV. 

The condensates of the operators in the vacuum are (uu)q = {dd)o = 
-(250 MeV) 3 and (^G 2 ) = (330 MeV) 4 . § In the QCD sum rule for the 
nucleon in the vacuum, ujq is determined to be 2.2 GeV by Borel stability 
analysis. 

u) + and w_ are determined by Borel stability analysis. Namely, we search 
for the values of cj+ and uj- for which the calculated strength a has the most 
stable plateau as a function of the Borel mass Mb both in the triplet and 
singlet channels. We find that the optimum choice is w+ ~ 1.3 GeV and 
w_ ~ 1.3 GeV. Figures 1 and 2 display how sensitive the Borel stability is 
to U+ and u— Figure 1 displays a vs. the Borel mass for tu + = 1.3 GeV 
and Li- = 1.2 GeV, 1.3 GeV and 1.4 GeV. Figure 2 displays the same for 
uj_ = 1.3 GeV and uj+ = 1.1 GeV, 1.3 GeV and 1.5 GeV. One sees that the 
Borel stability is much more sensitive to cj_ than to uj + . This is because the 
<jj- dependence of the coefficients of the dimension-three and dimension-four 
operators is much stronger than the lo+ dependence. 
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FIGURES 
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FIG. 1. The calculated strength a as a function of the Borel mass squared, Mg. The solid 
lines represent the spin triplet channel, and the dotted lines represent the spin singlet channel by 
fixing uj + to 1.3 GeV. The top, middle and bottom lines are for oj_ = 1.2 GeV, u;_ = 1.3 GeV and 
oj_ = 1.4 GeV, respectively. 
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FIG. 2. The calculated strength a as a function of the Borel mass squared, M\. The solid lines 
corresponds to the spin triplet channel, and the dotted lines corresponds to the spin singlet channel 
by fixing cj_ to 1.3 GeV. The top, middle and bottom lines are for uj + = 1.1 GeV, uj + = 1.3 GeV 
and lo + = 1.5 GeV, respectively. 

The a is determined by taking the maximum values for u> + = 1.3 GeV 
and w_ = 1.3 GeV as 

«L = 1-7 fm, 



We see that the spin-dependent part is rather smaller than the spin- 
independent part. This is due to the small matrix elements of the spin- 
dependent operators. From Figs. 1 and 2 we should expect errors of about 
30% due to the choice of uj + and w_ . The errors due to the uncertainties of the 
nucleon matrix elements of the dimension-three and dimension-four operators 
are about 10% and 20%, respectively. Combining all the errors together, the 
above results have errors of approximately 40%. 

Let us compare the above results with experimental facts. Experimentally, 
the scattering length and the effective range have been found to be 

a 3 pn = -5.39 fm, r 3 pn = 1.75 fm, 
a* = 23.7 fm, r\ n = 2.73 fm. 



In the spin-singlet channel, the scattering length is so large that the 
strength, is expected to be approximated well by 

"pn ~ r pn/ 2 = 1-37 fm, 
which is in rather good agreement with the calculated result. 
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In the spin-triplet channel, the scattering length is not small. However, 
it is not as large as in the spin-singlet channel. Therefore, it is not easy to 
estimate the strength c? pn from experimental observables. However, there is 
a loosely bound state, deuteron, in the spin-triplet channel, while there is an 
almost bound state in the spin-singlet channel. This implies that the interac- 
tion in the spin-triplet channel is stronger (but not very stronger) than that 
in the spin-singlet channel. This tendency is consistent with the calculated 
results. 

III. SUMMARY AND OUTLOOK 

In this paper we have studied spin-dependent nucleon-nucleon (NN) in- 
teractions in the QCD sum rules. The basic object of our study is the spin- 
dependent nucleon correlation function, whose matrix element is taken with 
respect to the one-nucleon state. 

The dispersion integral of the correlation function around the nucleon 
threshold has been investigated in detail. The integral is given by the sum of 
the threshold contribution, which is due to the second-order pole term pro- 
portional to the scattering length, the continuum contribution, and the bound 
state contribution. When the interaction is weak, the integral is dominated 
by the threshold contribution and is proportional to the scattering length. As 
the interaction becomes stronger, the continuum contribution also becomes 
important. When the interaction is just so strong as to form a bound state, 
both the threshold contribution and the continuum contribution diverge, but 
their sum is finite. Thus the sum of these two contributions is given by the 
same form as in the case of the weak interaction. However, in this case the 
scattering length is replaced by one half of the effective range. Based on 
this observation we have defined the NN interaction strength through the 
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dispersion integral around the nucleon threshold. 

In the OPE of the correlation function, new operators, such as qj^q, 
qi5C^uQ, have to be taken into account. These operators do not vanish when 
the matrix element is taken with respect to the spin-nonaveraged one-nucleon 
state. We have calculated the Wilson coefficients of such operators. 

The obtained sum rules relate the spin-dependent NN interaction 
strengths with the spin-dependent nucleon matrix elements of the quark- 
gluon composite operators. The spin-dependent nucleon matrix elements such 
as {qi^lbq)N and (q'jbcr^q}]^ are related with the spin-dependent structure 
functions of the nucleon, g\, gi and h\. We found that the interaction strength 
in the spin-singlet channel is weaker than in the spin-triplet channel, but that 
the spin-dependent part of the interaction strength is considerably smaller 
than the spin-independent part. Experimentally, it has been found that there 
is a loosely bound state, deuteron, in the spin-triplet channel, while there 
is an almost bound state in the spin-singlet channel, which implies that the 
interaction is slightly stronger in the spin-triplet channel than in the spin- 
singlet channel. This seems to be consistent with the sum rule result. In 
the spin-singlet channel, the scattering length is so large that the interaction 
strength can be estimated by using the observed effective range, while in the 
spin-triplet channel the absolute value of the interaction strength is difficult 
to obtain from observables. The empirical interaction strength thus obtained 
in the spin-singlet channel agrees rather well with the sum rule calculation. 

The method used in the present paper can be extended to other hadron- 
nucleon channels. In particular, it is straightforward to apply it to the 
hyperon-nucleon channels. The obtained sum rules would relate the hyperon- 
nucleon interaction strengths with the nucleon matrix elements of the quark- 
gluon operators which include strange quark operators in addition to up and 
down quark operators. Also, the present sum rules can be used in the oppo- 
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site way. Namely, if one knows detailed information on various spin-dependent 
hadron-nucleon interaction strengths, one can obtain spin-dependent matrix 
elements of quark-gluon operators with respect to the one-nucleon state. This 
provides us with information such as the spin content of the nucleon. 
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